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<N Abstract 
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Matrix model describing the anomalous dimensions of composite 
operators in M = 4 super Yang-Mills theory up to one-loop level is 
<^ | considered at finite temperature. We compute the thermal effective 

action for this model, which we define as the log of the partition 
function restricted to the states of given fixed length and spin. The 
Qh! result is obtained in the limit of high as well as low temperature. 
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1 Introduction 

The gauge/string correspondence has a long history starting from the 't Hooft 
paper [1], showing that the perturbative expansion of a gauge theory can be 
organized according to the geometrical genus expansion for the Feynman di- 
agrams. The idea behind this is that string theory may serve as an effective 
description for a strongly interacting gauge system. The physical picture mo- 
tivating such a description is given by condensation of the gauge field flux into 
tiny tubes leading to a linearly growing potential between quarks/ ant iquarks 
and string-like behavior. 

The development of string theory lead to the conjecture of AdS/CFT 
correspondence [2] (see the classical review on the subject [3]). According 
to this conjecture the string theory can be described in terms of gauge fields 
which are nothing else that the collective coordinates of D-branes, the non- 
perturbative objects on which the fundamental strings can end. Therefore, 
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the gauge field description of strings makes sense in the limit of strong string 
interactions. 

In this picture the rank N of the gauge group corresponds to the number 
of D-branes. In the 't Hooft limit of large N and small gauge coupling <?ym 
such that gyMy/N remains finite, the above duality relation is reduced to 
the correspondence between non-interacting IIB superstrings on AdSs x S 5 
background and M = 4 supersymmetric Yang-Mills theory (SYM) on the 
Minkowski space which is the conformal boundary to the anti-de Sitter space. 

An important feature of the AdS/CFT correspondence is that it is a 
true i.e. Ising-type duality which means that the strong coupling dynamics 
is mapped to a weak coupling one and viceversa (see [4] for a review of 
these type of dualities in lattice theories and spin systems). This property 
of AdS/CFT correspondence, beyond having a huge predictive force, also 
prevents from any direct proof of the correspondence itself, since for such a 
proof one should solve at least one of the models at strong coupling. 

Last years a considerable progress was registered in the study of both 
string and gauge theories in the way of approaching the AdS/CFT conjecture. 
Among this an important point is the discovery of the integrability (see 
e.g. [5-8]). 

It was a breakthrough to realize that there are certain limits in which 
both string theory and SYM are reachable in the framework of perturbation 
theory. On the string theory side these limits are special geometrical limits 
while on the gauge theory side they correspond to particular subclasses of 
SYM composite operators for which one can extend the applicability of the 
perturbation theory [9-11]. 

All this time the analysis was performed at infinite N, which in partic- 
ular allows one to use the integrability [7,8,12]. Since 1/N corrections are 
expected to break explicitly the integrability, during this development they 
got much less attention. 

In [13-18], however, it was shown that one can map the anomalous di- 
mension operator on M = 4 SYM into an interacting spin system taking into 
account the finite N rank of the gauge group. The finiteness of N, results 
in the chain splitting and joining interaction. This process mimics the string 
interaction^] 

There is, however, still another approach to the same problem. It con- 
sists in the interpretation of the operator of anomalous dimensions as the 
Hamiltonian of a gauged matrix model [19-21]. Spin chains in this picture 

2 To be precise, the above approach neglects the trace identities for U(iV) matrices. 
This means that only polynomial 1/N contribution is taken into account and not the 
exponential corrections. This is very similar to the Quantum Field Theory, where non- 
analytic corrections are systematically dropped in perturbation theory. 
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correspond to gauge invariant collective states. The problem with such a 
description is that there are not enough developed tools for the study of such 
type of matrix models independently of the spin chain results. 

One possibility is to describe the model in terms of a noncommutative 
gauge theory by expanding the matrix fluctuations about an appropriate 
classical solution [19]. This approach however, has its own drawbacks since 
it leaves a number of unanswered questions regarding the arbitrariness of 
such a description as well as regarding the spin chain/string interpretation. 
In particular, in the noncommutative description the matrix model of it is 
not very clear what is e.g. a spin chain. 

The present paper aims to fill, at least partially, the gaps in the under- 
standing of the matrix model for the anomalous dimensions of M = 4 SYM. 
Here we consider statistically large but finite values of N. This means that 
we do not exclude the string interactions from our description. According 
to BMN analysis [9], the string interaction is rated as L 2 /N, where L is the 
aggregate length of all spin chains. This interaction becomes important when 
one considers long enough operators. 

In our approach we find it convenient to put the matrix theory in a finite 
temperature background. In the Yang-Mills theory this temperature corre- 
sponds to the Euclidean compactified Liouville or scale direction. According 
to AdS/CFT correspondence this is a compact time cycle in the AdS space. 
As it was pointed in [22] , the thermalization of AdS space occurs due to the 
presence of black holes. 

The Af = 4 SYM theory is a conformal invariant theory, therefore its 
dilatation operator can be identified with the Hamiltonian dual to the ra- 
dial time. Since the dilatation operator and the SYM Hamilton are related 
through a similarity transformation, the spectrum and therefore, the ther- 
modynamical properties of the matrix model should be the same as those of 
the thermal Yang-Mills which were considered in [23-25] There is a differ- 
ence, however. In the contrast to the above cited works where the thermal 
Yang-Mills was considered on S 1 x S 3 , the radial time temperature should 
correspond rather to SYM compactified on S A . 

On the other hand, a similar partition function was considered in [29-31]. 
In particular, in [31], the thermal partition function of the gauged oscillator 
is derived by counting the trace states using the Polya Enumeration Theorem 
(PET). The problem of such approach is that, strictly speaking, it requires 
iV to be infinite. Because of this the analysis can not be extended beyond 
the Hagedorn temperature. In contrast to this approach we compute the 
path integral for the matrix model rather then just counting the states. The 

3 The SU(2) sector can be regarded as a particular limit of TV = 4 SYM [26-28]. 
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advantage of our approach is that it allows a conceptually simple (although 
technically evolved) extension to higher loops in the spirit of "conventional" 
perturbation theory. It appears that the description we obtain is in some 
sense complimentary to the one of the PET approach: for very long opera- 
tors we get a N 2 scaling of the extensive thermodynamical functions, which 
signals for the string bit phase, while for the operators below critical length 
the thermodynamical potentials scale as N°, which is compatible with N- 
independent description obtained in the PET approach. 

The plan of the paper is as follows. In the next section we review some 
basic properties of the matrix model: we rewrite it in the real form in terms 
of Hermitian matrices and pass to the "second order formalism" . The sec- 
ond order formalism is important for the comparison the matrix model under 
study with the DLCQ-inspired BMN matrix model [9]. The comparison re- 
veals that in our matrix model there are additional terms which seem unlikely 
to be canceled by the higher loop correction or to be generated from higher 
orders in the BMN matrix model. 

Next, we describe the canonical quantization and conserved charges of 
the model. In the third section we analyze the spin chains as collective 
states of the matrix model. In particular, we analyze the regimes at which 
the matrix model describes the gas of free spin chains. In section [3] we 
compute the free energy and the entropy (which we will like to call 'effective 
action') in the high temperature limit. During the computation we observe 
a couple of rather remarkable features of the model. One feature is that the 
matrix model substantially simplifies in the mode expansion which allows 
one to solve the model exactly at large N analytically in exponentials of the 
chemical potentials. 

We analyze the situation in both cases of small and large chemical po- 
tentials and find out that for large enough value of the potentials (above the 
Hagedorn one) all terms scaling as iV 2 and iV cancel out leaving us with the 
next leading contribution^) On the opposite side, when the chemical poten- 
tials are small the model behaves as a system of N 2 interacting particles. 
In section H] we consider small temperature limit. In this limit the typical 
behavior is one of a iV-particle system. Finally we give the string theory 
interpretation and discuss the results. 



4 An approach based on random walks allows one to find this contribution too [32] . 
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2 Matrix model 



2.1 Dilatation operator and anomalous dimensions of 

M = 4 SYM 

According to the AdS/CFT conjecture, the string states correspond to com- 
posite operators of the gauge theory. Furthermore, identification of the re- 
spective charges in the symmetry groups on both sides puts into correspon- 
dence the string energy levels to the eigenvalues of the dilatation operator of 
J\f = 4 SYM. Classically, the value of the dilatation operator is given by the 
composition of the dimensions of its elementary components. In the quantum 
theory this is corrected by the anomalous dimension contribution. It also 
appears that the dilatation mixes different operators, which form, however, 
invariant classes which do not mix. A particular example which we analyze 
in this paper is given by the SU(2)-sector of SYM, which consists of all gauge 
invariant polynomial operators built from two complex scalars, 

$ 1 = Z = i(0 5 + i0 6 ), (2.1) 
$ 2 = (t> = I(0 1 + i0 2 ), (2.2) 

where 0;, % — 1, 2, . . . , 6 are the scalars of the M = 4 SYM. 

The mixing matrix/dilatation operator in the context of the perturbation 
theory was considered in [8,12]. 

In particular, the one-loop contribution to the mixing matrix for operators 
in SU(2) sector was found to be [33], 

H = : tr[<T , $ b ] [$ a , $ 6 ] := -|m : tr[0, Z] [0, Z] :, (2.3) 

where the checked character corresponds to the differential operator, 

= 9&T> ( 2 - 4 ) 

and the colon ":" denotes such an ordering in which no derivative acts on 
the operators within the same group. 

2.2 The matrix action 

Consider the matrix model describing the anomalous dimensions of M = 
4 super Yang-Mills model. As it was found in [19] this matrix model is 
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described by the action 



S(%%A) = 

j dt ^tr^(^ a Vo^ a -Vo^ a )-tr^ a + ^tr[^ a ,^][^,^]^ 

(2.5) 

where^l \I/ a , a = 1, 2 are (non-Hermitian) N x N matrices, and Vo^ = do^ + 
[Aq, \&], with Hermitian matrix Aq, is the covariant time derivative. 

2.3 Real form 

To compare our model with other matrix models arising in the context of 
AdS/CFT correspondence it is convenient to pass from the variables \l/ and 
\l/ to the Hermitian ones X and Y representing the Hermitian and the anti- 
Hermitian parts of \I/ as follows, 

V a = j= 2 (P a + iX a ), % a = j-{P a -iX a ). (2.6) 
Substitution of ( 12. 6R into the action ( 12.51) yields, 

S(P,X,A) = 

j dt (tr \{P a {VoX a ) - (V P a )X a } - \{P 2 a + Xl) 

+ §^tr{[P a ,P b } 2 + [X a ,X b ] 2 + 2[P a ,X b } 2 }y (2.7) 

where we used the constraint [P a ,X a ] = and dropped off the total time 
derivative term (i/4)^(P a 2 + X 2 a ). The action (J23D can be written in a more 
economical way using a common variable X A of extended dimensionality 
A = 1,...,A, 

X A = (P a ,X b ). (2.8) 
In this case the action (12. 7p can be rewritten as 

S(X,A) = j dt (luj AB trX A V X B -±trX 2 A + §^tr[X A ,X B } 2 y (2.9) 
5 Note, however, the difference in the notations. 

6 In [19] these fields are denoted by X and X. Here, for convenience, we replace them 
by iff and iff, respectively. Also we introduce the term tr^vE' which describes the classical 
dimension. 
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where the anti-symmetric matrix ujab is given in terms of (X, Y) decompo- 
sition of Xa by the following block structure 

Let us note that the form (12.91) of the action is very close to the one of 
the "classical" Yang-Mills type matrix models. The only difference is the 
mass term and the first order kinetic term. 



2.4 "Second order" formalism 

The existend matrix models related to non-perturbative string dynamics are 
second order acrtion i.e. their kinetic terms have the form of velocity-squared. 
Our model is described by first order action i.e. with at most linear velocity 
dependence. On the other hand the matrix variables are non-Hermitian 
matrices, from which we can deduce that the action (12.51) the "first order 
form" of higher order action. Indeed, the example of ordinary harmonic 
oscillator with the classical action 

S = J dt (\i 2 -\x 2 ), (2.11) 

can be equivalently rewritten in the canonnical form 

S= J dt (px - \(p 2 + x 2 )) , (2.12) 

by the Legendre thransform with canonical momentum p = dC/dx. Further 
one can pass to the complex coordinate a which is the classical counterpart 
of the oscillator annihilation operator 

a = ^(p + ix), d=^(p-ix). (2.13) 

In the last coordinates the action takes the form 

S = J dt (| (ad — da) — da) + boundary terms. (2-14) 

If we replace a with the matrix valued ^ and take the trace we see that the 
quadratic part of the action ( 12. 5ft is nothing else that the harmonic matrix 
oscillator in the complex first order form, while ( 12. 7ft gives the ordinary first 
order formalism. Now, let us try to reconstruct the corresponding "second 
order" action. Classically, one can switch back to the Lagrangian formalism 
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just by solving equations of motion with respect to p in terms of x and x and 
substituting p in the action by this solution. 

The equations of motions resulting from the variation of P a in (12.71) reac|§ 

P a = X a - fp ([A, [A, P a ]] + [X b , [X b , P a }}) . (2.15) 

We can not solve this equation exactly, but we do not need the exact solution 
either. Recall that the action (12.51) is only the order g Y M approximation of 
the "exact" matrix model describing the anomalous dimensions to all orders 
in perturbation theory as well as nonperturbative effects. So, for us it suffices 
to solve the equation (I2.15P to only the first order in g\ M . This solution is 
given by, 

Pa = X a -&£ {[X b , [X h X a ]] + [X b , [X b , X a ]]) . (2.16) 

Then, the resulting "second order formalism" action takes the following 
form (up to the order Jym), 

S = 

J dt tr (|(X 2 - X 2 ) + §p ([X a ,X b ] 2 + 2[X a ,X b ] + [X a ,X b ] 2 )) . (2.17) 

The Gauss law constraint at the same time takes the following form, 

G = [X a , X a ] + ^ ([X b , [X b , X a ]] + [X b , [X b , X a ]]j . (2.18) 

Let us note the fact, that for slow varying matrices: X — > 0, the matrix 
action (12.171) is equivalent to the matrix model of [34]. This is not surprising 
since the respective matrix model reproduce the BMN spectrum, while our 
model reproduces the leading order in A of this spectrum. The surprising 
part is that beyond the ordinary commutator term our action contains also 
commutator of velocities. One can see that such terms can not be canceled 
by higher loop contribution. It would be interesting to see wether these terms 
can be removed by a proper redefinition of the fields. 

2.5 Quantization, gauge invariance 

To quantize the model given by the action (12.51) we have to impose first 
the gauge fixing condition. The most natural condition in our case is the 
temporal gauge A = 0. 

7 In quantum theory in addition to this one should also take care on modification of the 
measure (or of the scalar product). 
8 We use the gauge A = 0. 
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As it can be seen, the variation of the action with the respect to the gauge 
field A fails to produce a bona fide equations of motion, one has instead the 
Gauss law constraint, 

G = = \u AB [X A ,X B ] « 0. (2.19) 

From the simplectic structure of the classical action (12.51) one can extract 
the canonical Poisson bracket of the model. In components it reads: 

{(# a y", {^ h ) k 1 } = -i^V. (2.20) 

This can be written in the componentless form using the "quantum group" 
notations 

^ 2) } = -iP 12 , (2.21) 
where the N 2 x iV 2 -dimensional matrices \td and Pu are defined 

*W = tt <g> I, = I ® * a , (2.22) 

while P12 is the permutation operator 

P l2 a®b = b®a. (2.23) 

The canonical quantization consist in the promoting of the Poisson bracket 
in either form ( 12.2011 or ( 12.211) to the quantum commutator for which 

p a y , (%) k l ] = SX*. (2.24) 

The interpretation of the quantum algebra ( 12.241) is the following. The op- 
erator (^f a )^ creates a matrix element in the row i and column j having the 
polarization a while the respective component of destroys it. Let us 

note that the quantum commutator is denoted by the double braces in order 
to distinguish it from the ordinary matrix one which is related to different 
contractions of the matrix indices i,j, . . . but not related to the permutation 
of the operator valued matrix elements. 

In the complete analogy with the Harmonic oscillator the (extended) Fock 
space can be constructed by the action of the rising operators (^ a V on the 
oscillator vacuum state 

|c> = ' • • • |n> , (2.25) 

k 

where 

(* o V'|J2)=0, Va,z,j. (2.26) 
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Due to the constraint (12.191) the extended Fock space is too large. One 
has to restrict it to the gauge invariant subspace by imposing the condition, 

G|$)=0, (2.27) 

where G =: [$? a , \l> a ] : is the quantum version of the constraint f)2.19p . (Where 
the colon denotes an ordering prescription according to which all the com- 
ponents of \l/ stay always to the right of those of \l>.) 

As in the case of ordinary gauge fields the operator G corresponding to 
the constraint is the generator of gauge transformations: 

[G(u), M>J = -[it, [G(u), V a ] = [it, ttj, (2.28) 

where G(u) is the operator G smeared with an element of su(N) algebra u: 

G(u) = tr Gu = tr : [* a , a ] : u. (2.29) 

As a result of application the condition (12.271) the physical space is gen- 
erated by only those \C) which correspond to gauge invariant polynomials of 
\l/ a . In terms of equation (I2.25P this means that the quotients (C ,ai '" afc )} 1 1 '.'.'.^ 
should be linear combinations of the products of delta symbols <V . The most 
generic such monomial of order L is given by a permutation of L elements 

7 e S L , 

(C7 Bl -- 0fc )£;;** = C 01 "- *' 7 ^ . . . 5 jk 1 ^, (2.30) 
where 7^ = ^(k) is the permutation of the A;-th element. 

2.6 Conserved charges 

In the model (I2.5P there is a number of quantities which are conserved. 
Among these let us consider the quadratic ones 

L a6 =trf a * 6 , L[ b = L ba , L ab = 0. (2.31) 

They can be split into Hermitian operator L and operator valued vector S, 

L = L aa , S = <r a bL ab , (2.32) 

where a is the vector of Pauli matrices. As we will see below these operators 
correspond to the total length and total spin operators of the spin chains. 
Both L and 5* generate a representation of the algebra U(2) which can be split 
into irreducible components described by positive integer value of L (U(l) 
spin) and the spin s irreducible representation of the SU(2) component. 

In particular, due to the normal ordering of the operators all above 
charges vanish on the vacuum state The value of the operator L cor- 
responds to the total number of the excited oscillator modes, while e.g. 
5*3 = \ tr ^i^/i — ^2^2 corresponds to the difference in the number of excited 
modes for the first and the second oscillator. 
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2.7 Matrices as spin chain gas 

As we discussed earlier, the quantum Hamiltonian of the matrix model can 
be mapped to the Hamiltonian of a spin system with a chaining interaction 
[13,14] (see also [20,21]). The limit iV -> oo with A = # 2 iV-fixed reduces 
that model to the integrable Heisenberg XXXi/ 2 model. In this picture spin 
chains correspond to one-trace operators. As soon as N is finite the subspace 
of such operators is not invariant with respect to dynamics and, therefore, in 
a gauge invariant configurations one is forced to consider a mixture of matrix 
states with all possible number of traces. 

Naively, in N — > oo limit the non-planar interaction vanishes and one 
ends up with an "ideal gas" of spin chains, where each chain is conserved by 
the dynamics. 

This simple picture, which we have for the spin chain/matrix model cor- 
respondence in the quantum theory, breaks down as soon as we consider the 
(semi) classical limit of either spin chains or matrix theory. As it turns out, 
the semiclassical regimes for two descriptions appear at different if not con- 
tradictory conditions. Indeed, the matrix model when the expectation values 
for the occupation numbers of most modes are larger than one. This gives 
us the condition, 

L = tr$A ~ iV 2 . (2.33) 

On the other hand, as we know from the BMN evaluation, which should hold 
also in our case, the non-planar interaction rate is given by a factor propor- 
tional to L 2 /N. Then, the regime (12. 33ft for the spin description means that 
the non-planar interactions are not only strong but overwhelmingly dominat- 
ing the dynamics. In this case one can not even speak about the spin chains. 
In contrast, for L <C y/N, when one can neglect the non-planar interac- 
tion and, hence, the spin chains are stable the matrix model is in essentially 
quantum regime. 

A gauge invariant matrix state generally will contain a mixture (gas, 
condensate, etc.) of spin chains of different lengths. The dynamics of this 
"spin soup" is such that the very short chain will tend to join bigger ones 
while too long chains will tend to split into smaller ones. Therefore, an 
equilibrium state should be statistically dominated by chains of a particular 
length. 

As one can see, the statistical description of such a complicate system 
appears very naturally as it can catch the properties of the dynamics of both 
spins and matrices. Moreover, moving from one regime to another we can 
extract the relevant information like which sort of behavior dominates at 
different regimes. We do this in the next sections. 
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3 Spin chain gas at high temperature 



Let us consider the closed sector of matrix model given by states of the length 
L, and let us further restrict ourself to the states of definite value of spin S. 
The partition function restricted to such a subspace of the Hilbert space is 
given by, 

Z{L,S) =trn(L,5 ? )e- /3H( *'* ) , (3.1) 

where I1(L, S) is the projector to the (L, S)-subspace and H(^/, \P) is the 
Hamiltonian operator of ( 12. 5ft . The projection to given values of L and S 
can be realized in the following way: 

Z(L,S) = J dfidxe^ L+ ^Z(fi,x), (3.2) 

where Z(fi, x) is the (grand canonical) partition function: 

Z(n,x) = tre- ff ^*'*), (3.3) 

where the chemical potentials fi and x enter into the Hamiltonian H^(^f, 
in the following way, 

t r r-i* a [A*a]+ A*a+5^«^«6*6-^^[*«,*6][*^* 6 ^ ■ ( 3 - 4 ) 

The original Hamiltonian is recovered at fi — (3 and x — 0. 

To find the partial partition function ( 13. 2p let us first compute the grand 
canonical partition function Z(fi,x). 

An important issue in this computation is one of gauge symmetry and the 
gauge fixing. There are several possibilities with this. One possible approach 
is given by imposing the diagonal gauge fixing to one of the components, 
say the real part of As is claimed in [34] this leads to the description 
of the model in terms of a fermionic liquid. In the present case we use 
a different approach. To discuss the problem of gauge fixing it is useful 
to switch for a while to the path integral interpretation of the partition 
function. In this interpretation the inverse time (3 plays the role of the size 
of compact Euclidean time. Therefore, the path integral formulation for 
the partition function implies periodic boundary conditions for the matrices: 
X(t + (3)=X(t). 

Let us now turn to the gauge field. It is clear, that because of periodic 
time we can not impose the temporal gauge: A = 0. Indeed, out of the gauge 
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field A(t) one can form a gauge covariant object called Polyakov loopo 

W[A] = Texp (j) drA(T), (3.5) 

whose eigenvalues are gauge invariant and, therefore, can not be canceled by 
a gauge transformation. Hence, the admissible gauge we can impose is to fix 
A to be at most constant and diagonal. 

In the operator formalism one can still overcome the above restrictions 
and impose the temporal gauge A = 0, but one has supply it with the Gauss 
law constraint on the matrices. As the Gauss law constraint is a conserved 
quantity, it suffices to impose it once at a preferred time instant. This is 
equivalent to the restriction of the gauge field A in the Hamiltonian ( 13. 4j) to 
be constant and diagonal. 

The Gauss law constraint act on the quantum states as generator of gauge 
transformations. Hence, vanishing of the constraint is equivalent to the gauge 
invariance of the respective state. As we discuss in the Appendix (Bj in 
the case of compact Lie gauge group, projecting to the subspace of gauge 
invariant states can be implemented by the group integration. The advantage 
of this approach is that the compactness of the gauge group is taken into 
account as well. 

At large temperature and fixed value of the length operator L the commu- 
tator term of the Hamiltonian ( 13. 4p is suppressed by a factor /3gy M , therefore 
we can treat it as a perturbation over the quadratic part of the Hamiltonian. 

In what follows we will use the holomorphic representation which cor- 
responds to the quantum oscillator described by the quadratic part of the 
Hamiltonian ( 13. 4ft . For the convenience of the reader we give a brief summary 
of the holomorphic representation in the Appendices [A] and [B] 



3.1 Quadratic part 

As we discussed above, we can consider the commutator term to be a small 
perturbation over the quadratic part of the Hamiltonian, such that we can 
apply the perturbation theory. Let us first compute the bare partition func- 
tion which is the trace over gauge invariant subspace, 

Z (P)=tr gA .e- H °^, (3.6) 

where ^ is the quadratic part of the Hamiltonian (13. 4p . which consists 
of chemical potentials only: 

H 0> ^(V, *) = tr (n^a^a + \^> a {x ■ a) ab V b ) . (3.7) 

9 This equation defines W[A] up to a gauge transformation related to the choice of 
initial point. 
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The Gaussian integral in (13. 7p is convergent provided \fx^ < /x. For lager 
values of x it diverges because the quadratic form in (13. 7p becomes indefinite. 
Physically, however, more interesting situations are when spin is suppressed 
stronger than the length operator. In principle this can be done by a "Wick 
rotation" to imaginary spin chemical potential, x \— > ix. 

Restriction of the trace to the gauge invariant subspace can implemented, 
as discussed in Appendix |Bj by averaging over the global gauge group, 

J dU J exp (- ti y aU- 1 ^ a U + fe^**/ 2 *) , (3.8) 

where the exponent in the last term is understood as a 2 x 2 matrix exponent. 
Let us note that, in contrast to ordinary gauge fixing, averaging over the 
gauge group keeps the gauge invariance explicit i.e. the integral is invariant 
with respect to global gauge transformations: 

U -> V~ X UV, -> V' 1 ^, f -»■ V~ l W. (3.9) 

In the basis of eigenvectors |±) of the matrix x ■ a, 

x • a |±) = ±x |±) , x = Vx^, (3.10) 

the matrix e -^~ x - a / 2 becomes diagonal with values e~ M± , where 

x 

fJL±=fJL±-. (3.11) 

The symmetry which is left allows us fixing the matrix U to be diagonally 
Also changing the variable in the [/-integral to the diagonal values U = 
diag{e ie "} produces square of the Vandermonde determinant, 

\A(9)\ 2 = 11 2(1 - cos# mn ), e mn = 9 m - 6 n . (3.12) 

We also used the fact that the integral does not depend on the U(l) factor. 
Hence one can fix the sum of the eigenvalues ('center of mass') to vanish f^l 

J29 n = mod 2tt. (3.13) 



10 This is in the direct relation to the fact that the gauge field can also be fixed to be 
constant diagonal. 



li 



In what follows we just systematically drop the contributions proportional to ^2 On 
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As a result of above transformations the gauge invariant trace takes the 
following 'normal form', 



Z {v,x) = J Y[d9 n \A 



I \ 



x / ~^ ex P 



(3.14) 



■ m,n i 
\ ± / 



where the normal modes h^Jfi, x, 6), are given by, 

= e Wm " - . (3.15) 

As one can see, the integral over the matrices \I/ and * is Gaussian and 
can be easily taken. The quadratic partition function then is reduced to the 
integral over #'s only, 



m,n 



2 -fJV(JV+l) e 7VV 



[sinh(/z+/2) sinh(//_/2)] 



N 



U d9 « 



TT 1 - COS fl mn . 

(cosh /i + - cos 6 mn ) (cosh /i_ - cos 6 mn ) ' 

Using ( 13.21) . one can pass to the variables L and S 1 . As soon as a large 
number of degrees of freedom is concerned this integration can be done by the 
saddle-point approximation, i.e. by replacing /x and x with their solutions in 
the presence of the sources L and S to the equations of motion corresponding 
to the effective action resulting from (13.161) . The saddle point conditions 
coincide with the Legendre equations. Hence the Laplace transform in this 
limit is reduced to the Legendre transform of the free energy: 

S cS (L,S) = Ln(L, S) +iS ■ x(L, S) - F(n(L, S),x(L, S)) (3.17b) 

where functions /i(L, S) and x(L, S) in the second line are found as solutions 
to the first line equations ( 13.17aj) . The free energy F is defined as the log of 
the partition function, 

F(/i, x) — — lnZ (//, x). (3.18) 
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Thus, to obtain the effective action we have to integrate over the gauge field 
eigenvalues 9 n . In some cases this can be also done by the saddle-point 
approximation, i.e. replace the integral by the value of the integrand at its 
maxima. 



3.2 The free energy 

Consider the saddle point approach to the integral ( 13.16p . Finding the max- 
imal value of the integrand in the partition function (I3.16P is reduced to the 
problem for finding extrema with respect to variation of 9 n for the following 
action, 

F{9; n, x) = -N 2 fi + N[ln sinh(//+/2) + In sinh(/i_/2)] 
+§ 2^ ( _ l n (l _ cos9 mn ) + ln(coshyU + - cos6> mn ) + ln(cosh/i_ - cos9 mn )) , 

m,n 

(3.19) 

where we dropped the constant terms, which can be absorbed by a redefini- 
tion of the measure. 

The equations for the Legendre transform (13.171) then read, 

L = -N 2 + iJV[coth(//+/2) + coth(/x_/2)] 

+ \ £ ( , Sinhfl+ a + h Sinh ^ a ) ( 3 - 20a ) 
w V cosn /•*+ - cos C cosh /i_ - COS 9 mn J 

25= iJV[coth(/! + /2) - coth(/<_/2)] 

+ I V ( Sinh "- ) -, (3.20b) 

ln~?i \coshfx + - cos6> mn cosh/i_ - cos 9 mn J x 



E| sin 6> nm ^ sin 6> nm sin 
\ 1 — ct r< () r< wli ii , — i ■! >- 



cos 6* nm cosh /i + — cos 9 nm cosh //_ — cos 9 m 
= 0, (3.20c) 

where the summation in the last equation should run through m different 
from n. As some 9 m approaches the 9 n , the first term in (I3.20cl) diverges 
as ~ l/(9 n — 9 m ). In what follows we assume that this term is regularized 
by replacing 1 by by a term coshe, where e afterwards will be sent to zero. 
Having this regularization in mind we can extend the last sum to all m. 

Let us first discuss the equations ( 13.201) from the qualitative point of view. 
The last equation, (I3.20cj) describes an equilibrium static configuration of 
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the system consisting of N particles on a circle with the following pairwise 
interaction potential, 



ip = — | ln(l — cos A) + | ln(cosh fi + — cos A) + | ln(cosh \i_ — cos A). (3.21) 

As one can see, when \x± 7^ the interaction potential consists of one 
repulsive and two attracting terms. Normally, at small distances the repulsive 
term dominates while at large separations the interaction becomes attractive. 
Now several variants are possible. First, suppose the scale at which the 
attractive force starts to win over the repulsive one is larger than 2tt. In this 
case the repulsive dominant interaction will force the eigenvalues to arrange 
uniformly over the circle. This is a configuration with unbroken subgroup of 
the translation symmetry for 9: 6 n *—>■ 9 n + rj. This happens when n± is large 
enough. When n± decrease the scale at which interaction becomes attractive 
decreases as well and at some point it becomes less than the circle's length 
27r. At this stage for eigenvalues it becomes more "convenient", from the 
point of view of "total energy", to condense around some point (s) rather 
than be uniformly spread over the circle. In particular, when at least one 
of /a± vanishes the repulsion term is canceled completely and all eigenvalues 
tend to collapse to a single point. The last is the extremal case of breaking 
of the translational symmetry of the eigenvalue distribution. 

3.2.1 Expansion in powers of e~ M± 

Let us pass to a quantitative analysis rather a qualitative description. The 
problem dramatically simplifies if one looks for solutions which are analytic 
in e~ fl± . This appears possible due to the fact that the expansion in terms of 
powers of e _At± of the 'equation of motion' (I3.20p is remarkably simple and 
takes the following form, 



As one can see, this expansion organizes itself into something similar to a 
Fourier series. 

Applied to the free energy (13.191) the expansion yields, 



00 





F(6; /i, x) = N[\n sinh(// + /2) + In sinh(ji_/2) - ji] 



00 1 

+ V-(l-e-^+ -e- 



) ^ cos(u9 mn ). (3.23) 



m,n 
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In (I3.22p and (I3.23P one can observe that the terms scaling as N 2 were can- 
celed away by the zeroth term in the expansion. To obtain (13.221) and (I3.23P 
we used the following expansion, 



In (cosh p — cos6) = p — In 2 — cos[n6], (3.24) 



n=l 



as well as its derivative. 

Let us continue with the succession of cancelations. Observe the following 
properties of trigonometric sums, 

2^ sin uo9 mn = sin u9 n cos u9 m — cos u>9 n sin u9 m , (3.25a) 

m m m 

cosuo9 mn = I ^^cosu;^ J + I ^^sinct;^ J (3.25b) 



m,n 



where, 



for io = 1, 2, 3, . . . 



~ ( cos2 u ® m + sin2 u ® m ) = + ~ N ' 



^cos(u;6> n ), p u = sin(a;0 n ) (3.26) 



It is natural to assume that the equilibrium ^-distribution possesses at 
least one reflection symmetry such that one can choose a point with respect 
to which the distribution of #'s is even. In this case the sums over sines p^ 
vanish and we are left with the following equations 

oo 

(- 1 + 2e~^cosh (^)) wi(uOn)Pu, = 0, (3.27) 

w=l 

as well as with 

oo -. 

F{9- /*,5)=£-(l- 2e ~^ cosh (t)) ^ ■ (3 ' 28) 

U)=l 

A glance at (13.271) and (I3.28P reveals that the free terms scaling linearly in 
N are canceled if we extend summation to all possible pairs of eigenvalue 
numbers. 

It is interesting to note that equations of motion do not depend on the 
constant part of eigenvalue distribution p . Hence, an obvious solution to 
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the equation (13.271) is given by the uniform distribution of the eigenvalues: 
9 n = 2nn/N, n = 0, 1, 2 . . . , N - 1, or p(A) = N/2ix = constant. 

Let us find the other solutions. In the limit of large number of eigenvalues 
their distribution can be described by a continuous density p{6) such that, 



d9p(6)=N, p(6)>0, (3.29) 
and a sum over eigenvalues can be replaced by an integral as follows, 

£/(*»)= <fd0p(0)f(0), (3.30) 

n 

for some function f{9) which is smooth enough. In particular, p w in this 
limit become coefficients of Fourier series, 



Pu 



/oo _, 
dA p(X) cos(c<jA), p(A) = po + — coscjA. (3-31) 



Consider now the equation (I3.27P . The equality in (13.271) holds when each 
term in the sum vanishes separately i.e., 

p w (-1 + 2e-^cosh (jy)) = 0, w = 1,2,3.... (3.32) 

Eq. (13.321) in particular means that the Fourier expansion of p(9) can have 
only such modes (p u ^ 0) for which the expression inside the parenthesis of 
(I3.32p vanishes by itself, 

- 1 + 2e~^ cosh (if) =0, cu e Z+. (3.33) 

So, let us analyze the possible solutions to (I3.33p . Basically the equation 
(I3.33P admits at most one solution when p > | . Even in this case, due to 
the discrete nature of u> such a solution is possible only for particular pairs 
(/i, x). However, when p and x become small u can be treated as a continuous 
variable, therefore the set of pairs (p, x) for which one can find an integer 
solution becomes dense. 

In general, for each value of u there is a one-parameter family of (p^, xj) 
satisfying (I3.33p . 

Now, let us turn to the "free energy" F(p,x). Due to the factor (1 — 
2e _£JM cosh multiplied to each mode p u the "on-shell" value of F(9; p,x) 
vanishes. Physically this can be related to the following fact. For simplic- 
ity consider the well-studied case of the gauged one-matrix oscillator. This 
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system, unlike the multi-matrix oscillator, can be solved by imposing the di- 
agonal gauge to the matrix field. The Faddeev-Popov determinant, or better 
to say the Jacobian, which arises when passing to the description in terms 
of eigenvalues is again a Vandermonde determinant. Upon quantization the 
determinant can be absorbed into the wave function by rescaling it by the 
square root of the Vandermonde determinant. Such a rescaling makes the 
wave function odd with respect to permutation of the eigenvalues of the ma- 
trix field. Thus, the diagonal matrix oscillators behave like fermionic ones. 
In particular, they have to obey the Pauli exclusion principle: no two diag- 
onal oscillators could be in the same state. One of the effects of this is that 
the states with L < N are banned. In terms of the chemical potential this 
means that for \x larger than a critical value \i c the free energy should become 
trivial i.e. independent of \l. Rather than disappointing this fact seems to 
be compatible with the PET calculation of the partition function [31], which 
establishes that the free energy in this regime is finite (and N independent) 
as N goes to infinity. 

In spite that it appears that our approximation was not sensitive enough 
to obtain some nonzero thermodynamical quantities we still can try to ex- 
tract some information from the system by observing that writing down the 
Legendre transformation equations (13.201) before solving equations of mo- 
tions for #'s lead to a non-trivial result for L and S. This corresponds to the 
interchange of integration in #'s with one of \x and x. 

In the case of a single nontrivial mode ou, satisfying (I3.33p . the Legendre 
equations (13.201) are reduced to, 



where we used (13. 33ft to eliminate the factor 2e~ UJ/1 cosh(a;x/2) in the first 
equation and 2e _aJM in the second one. Equations (13.341) together with the 
mode equation (I3.33P seem to be enough to eliminate \l and x in favor of L 
and S for given uj: 




L 



(3.34a) 
(3.34b) 




(3.35b) 



(3.35a) 



As the "free energy" F vanishes on-shell the effective action S e s(L, S) is 
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given by, 



Suj = Lfi(L, S) + S ■ x(L, S) = 

+ 2Stanh~ 1 ( 1 I . (3.36) 





Next, according to the saddle point approximation, we have to sum over 
the contributions of all saddle points i.e. to sum over u>, 



e Seff 



J2 e ^- ( 3 - 37 ) 



Some comments are in order. The sum (13.371) diverges, but let us recall 
that the validity of our approach was given by the approximation of smooth 
eigenvalue distribution. This is definitely not true for the frequencies u > N, 
since they involve wavelengths shorter than the average distance between 
the eigenvalues 9 n . So, the applicability of the approach is restricted to still 
relatively large values of \i and x. Under relatively large we understand much 
bigger than ~ 1/JV. Also the equation (13.271) implies (13.321) only in the limit 
fj,N ^> 1. Indeed, for finite N the sines/cosines of u > iV can be expressed 
in terms of sines/cosines of smaller arguments. Therefore, (13.321) will get 
corrections from the higher modes: u > N. For fi ^> 1/N this contribution 
is exponentially suppressed and we can neglect the modes higher than N, 
otherwise they should be taken into account. As a conclusion the sum in 



(13.371) should be understood as regularized by restricting it to lu < N. 

It is clear, however, that as soon as large numbers are concerned, ^> 1 
and we have, 

e Si > e 5 u = ^j-^ Vw > 1. (3.38) 

Hence, for large L all terms except the first one can be neglected and we have 
SW = Si- 

As one can see from (13.361) . there is no iV and iV 2 scaling left in the 
effective action. It would be also interesting to check the effective action 
(I3.36P against the PET result. 



3.2.2 Solution for small [i and x 

Small chemical potentials correspond to larger expectation values for respec- 
tive thermodynamical quantities, since the suppression of configurations with 
large charges is weaker. Therefore considering small values of \x and x will 
favor contributions with large L and, respectively, S. 
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In the approach of previous subsection, however, this limit is singular 
since the expansion in powers of e _/1± diverges for small fi±. 

Now, since fi± are small we expect that the eigenvalues 9 n should condense 
around some arbitrary point breaking the translational symmetry. Let us 
pass to the quantitative analysis of this well. 

For small values of fi± the trigonometric and hyperbolic cosines can be 
replaced by the first terms of their Taylor expansion, 1 ± x 2 /2 + . . . . In this 
limit the action assumes the following form, 

S(6;v,x) = -iV 2 (/i-ln2) +iV[ln(/i + //_) -2 In 2] 

+ IZ)(- ln ^« + ln ^+ + d ^n) + M/*- + O) • (3-39) 

m,n 

The "action" (13.391) leads to the following "equations of motion" , 

mj L n \ °nm p+ ' u nm h 1 - ' °nm J 

Let us note, that when at least one of /i± vanishes, the potential becomes 
purely attractive and all eigenvalues collapse to a single point. In general, 
potential becomes attractive outside the region 9 mn < ^//JZjT/IZ - , so the eigen- 
values should condense to a region of this typical size. 

The problem of finding the eigenvalue configuration satisfying (13.401) is 
equivalent to finding an equilibrium distribution of particles on a line with a 
Van- der-Waals- like pair interaction given by the potential 

ip(X) = | (- In A 2 + ln(/4 + A 2 ) + ln(/i 2 + A 2 )) . (3.41) 

Again, in the limit of large N, the the eigenvalues become dense and one 
can replace the equation (13.401) by a continuous integral one: 

dvp(v) \~~\~~ + 2 *7, V , 2 + 2 ^TT v ) = °> (3- 42 ) 
V A -rj fi 2 + + {X-r]) 2 /ii + {X-r]) 2 J 

where the equation should hold for A G supp p(A). 

We approach this problem in the Appendix [D] The resulting 'energy 
function ( ID. 101) generates following Legendre equation^ 2 • 

L = 2iv 2 - p: 1 ) c d£ 



2S = 2N 2 (^ + ^) e| 



(3.43a) 
(3.43b) 



12 Note that here and on we replaced the "original" length L by a "renormalized" one 
L-N 2 . 



23 



where the function 8 is given by (see (ID.10I) in the Appendix [D]) 
= ^arctan^" 1 + ^ _1 arctan^, £ = 




~" (:Ui) 



In terms of L± = L ± 2S Legendre equations (13.43al) take the form, 

L± = ±4iV 2 /ii 1 ^ , (0- (3-45) 
Dividing the equation for "+" by the one for "— " , we obtain, 



l. 



(3.46) 



A remark is in order. Physically meaningful region is one where < 
L < L < L + . This corresponds to such a regime where chemical potentials 
statistically suppress the spin stronger then the length L. On the other 
hand, we know that the matrix integral diverges at fi < x/2. So, to reach 
the physically relevant region we are forced to make an analytic continuation 
to complex x and £ 2 < 0. 

Eq. (ET35j) together with (ET3S1) yield, 



// ± = ±4N 2 Lg 1 C£'(0, 
while in terms of \x and x this solution takes the form 

/ i = 2iv 2 (L; 1 -L: 1 )^'(0, 

x = 2iV 2 (LT 1 - LZ^S'it) 



(3.47) 



(3.48) 
(3.49) 



Plugging this solution into the Legendre transform equation (13.17bl) . we 
obtain, 



S e s{L, S) = |(L+// H 



AN 



where L± are given by, 



L_//_) - AN 2 S(0 = -4iV 2 ^(0 
L 




1 + Jt 



+ 



hi 



1 + JP- 



1- ^ ^ 



L±2S. 



(3.50) 



(3.51) 



Let us recall that e 5cff( - L ' 5 ^ gives the number of composite SYM operators 
of length L and spin S. It is instructive to compare eq. (13.501) to the log of 
the number of states of given spin in a system of L completely random spins. 
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3.3 One-loop contribution 

So far we considered the purely oscillator part of the matrix model. As we dis- 
cussed above in the limit of high temperature and small coupling the quartic 
commutator term can be considered to be a small perturbation. Therefore, 
in this limit one can find the correction to the oscillator effective action by 
just evaluating the average of the quartic term in the oscillator dominated 
background. This amounts to neglecting among others the back reaction of 
the background to the presence of the perturbation. At the one loop level it 
is a justified approximation, since the back reaction comes at the order g Y M 
and are of the same order as two-loop contribution from SYmJ^I 

In the previous section we considered the quadratic part (13. 7p of the 
matrix model Hamiltonian (13. 4|) . 

H = H + V = 

da 2 

tr(A* a +|5 -Vat A) -t^tr[V a , (3.52) 

At high temperatures the coupling (3g\ M is small so we can expand the 
partition function as follows 

Z = e" 5cff = tre- ( ^ 0+y) 

= tre~ Ho - tre~ Ho V + ■■■ = Z (l - (V) + . . . ) 

we-^-Mo, (3.53) 

where the mean value (V)o is defined in a standard way 

(V) = Zo 1 tTe- Ho V. (3.54) 

Thus, computation of non-zero coupling correction at the leading order re- 
sides in the calculation of the mean ( 13. 54ft of the following operator 

Bn 2 

V =17TT tr[* a ,* 6 p a ,tt 6 ]. (3.55) 
To do it let us use again the anti-holomorphic representation (see Appendix 

ED. 

Since the operator (I3.55P is a normal one its anti-holomorphic kernel is 
given by Eq. (1A.9[) . 

K V (*,Q) = MMe tr ^tr[<S Q , $ 6 ][$ a ,$ 6 ]. (3.56) 

13 Taking into account corrections from higher orders could be similar to one-loop renor- 
malization group improvement. 
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Since both the bare quadratic Hamiltonian and the perturbation are 
gauge invariant it suffices to insert the projection to the gauge invariant 
subspace only once. Therefore, the trace (13.54p is given by 

Z (V) = 

Mm j deA (6) j ^ ex P L y: + htt&d*£)\ x 

\ mn / 

tr[$ + e^+,$_e^-] [$+,$_], (3.57) 

where hmn was defined in (13.151) of subsection 13.11 

Making the substitution: hmn&mn — > $mi the integral (I3.57P transforms 
into 

(V)o = ^ [ ^?e-^e- 2 -trP + //,W),($_//,(-))][$ + ,$_], (3.58) 

07T Z J 7T JV 

where (&±/h ( - ±S) ) is a matrix given by the elements 

($±/h (±) ) mn = $tJh£l (3.59) 

and we used Zq = l/fl hS + >hS~> to cancel Zq in the l.h.s. 

This matrix integral can be computed taking into account that, 



rai a 2 bi b 2 — f d^dffi x 2 6i b2 \ -tr#*_ 
■ I felilfe222"iiWiTO2n2 / 7I -A r2 ^ kih k 2 h mini ^m 2 n 2 ) K 

S ai 1 8 a2 2 5 nikl 5 mi i 1 5 n2 k 2 5 m2 i 2 + 5 aib2 5 a2bl 5 ni k 2 5 mi i 2 5 n2 ki5 m2 i 1 . (3.60) 
Now, plugging (I3.60p into (I3.58P we get, 



/^5 , YM„-2« I 1 / ' 

n i 

P9ym ( N 



( V )o Jjfe ^ I h (+) h H 5 £ I h (+) h (~) ^ 

\ n Unn Unn knm \Hmnll> nk Hmnll> nk 



knm 



Aix 2 V(e M + - l)(e^- - 1) 
cos 6> mfc - e" M+ cos nfe - e _/i - cos 6 mn + e" 2/i 



(cosh/i + — cos 6 mn ) (cosh — cos# 



nfc ; 



<V)„+<V%. (3.61) 



The term in the first line of the r.h.s. of (13.581) can be readily evaluated 
to be, 

~ (3g YM N e-"+-"- 

(T/)o(/i ' " } = 4^ (l-e-. + )(l-e-^r (3 ' 62) 
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or in terms of the level uj solution to the Legendre equations (see (I3.35P ): 

the same quantity can be written as, 

<V%(L, 3) = 0- i ^ i -. (3.64) 



(i3s)--i)((iSs)--i) 



The a; = 1 term is a constant, 



(V)o = ^- (3-65) 

It is interesting to note that for large uj the mean value (I3.64p behaves 
like, 



(V)o = 

f3X \ uj 2 uj 



+ , , 2L ^ + •■■ • (3-66) 



The overall negative sign in front of uj 2 term makes the terms with large uj 
to be even more strongly suppressed. 

3.3.1 Expansion in powers of e~ M± 

Let us consider the second line of (I3.58f) and expand it in powers of e~ M± the 
same way as we did in the previous section. Remarkably, again the expansion 
is not only doable, but also appears in a simple form, 

(V)' = (V)o -{V% = 



YM 



|t2 ^e-(- +1 ^-^' +1 ^-^cosM mn + ^ fcn ), (3.67) 

u>,u)' rank 

where uj,uj' = 1,2, .. . enumerate the modes of the eigenvalue density dis- 
tribution. 

As a warmup it is not difficult to see that for homogeneous distribution 
of eigenvalues 6 n the average of each such term vanishes i.e., 

(^)o = (V%. (3.68) 
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Now let us turn to the situation where the eigenvalue distribution is 
'perturbed' by a non-constant mode u. 

A straightforward way to compute the expectation value (V)o would be 
to plug the eigenvalue distribution into ( 13. 67ft and get the answer. In the case 
of analytic solution we have the possibility of an indirect way of evaluation 
of (I3.67P even without knowledge of the details relating to the density p(A). 

Indeed, let us expand the cosine in (I3.67P and use the even character of 
the eigenvalue distribution (see the Appendix O) to get, 



\ V /0- 4?r 2 e Z^ e Pu,Pw>Pu,+w> ~ 4?i . 2 e 



\ 

, (3.69) 



\ & J 

where we separated explicitly the zero mode po — N/2tt of the eigenvalue 
distribution density. Using the property (I3.33P of the modes of the density 
one can readily evaluate the first term of the last equality of (13.691) to be, 

(V)' = -M^ £ (e"-+ + e"--) p ~pl = 



P9ymN -2 M >p ~2 _ @ X -2m t 
87r 3 2^P- ( 2vr )3 e 

p\ r r (is 



(2tt) 3 I V L 



2 \ uj 



L 1 - -r » ( 3 - 7 °) 



or for w = 1 it becomes, 



(V)^ =1 = -^(L>-4S>)/L. (3.71) 
By contrast, for uj — > oo, this is, 

As about the second term of ( 13.691) . it should vanish, since there is only 
one non-trivial mode uj. 
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3.3.2 Solution for small ji and x 

Consider now the case of small /i and x. In this case the only available 
approach is to plug the solution for the eigenvalue density ( ID. 91) into ( 13. 61 f) 
to evaluate (^)o- 

The term (^)o in this approach is of the order N°, and, therefore is small 
with respect to the terms {V)' . 

Let us turn to the evaluation of (V)' Q . The leading term in /i ± and 8 mn 
of the second line of (13.61 j) reads, 



(V)' 



09ym 
4tt 2 



Mm 

47T 2 



dAp(A) I dr]p(rj) j dv p(u) 

V+V- ~ (v~ A) (A - z/) 



x 



(^ + (^-A) 2 )((^- + (^-A) 2 ) 



(3.73) 



In the incompressible liquid approximation of the Appendix O, the inte- 
gral in ( 13.731) becomes, 



(vy = ^ 



(3X 



dXx 



2(2tt) 2 A j_ a 

A 2 -(A + i/i + ) 2 



In 



A 2 - (A - i/i , 



In 



A 2 - (A + i/x. 
A 2 - (A - i//_ 



//4+^A + A) 2 
1 l X + (A-A) 2 



In 



/i 2 +(A + A) 2 
/i 2 _ + (A - A) 2 



(3.74) 



where we use the shortcut notation A = ^fjn^TZ. Performing a substitution 
A — > AA in (13.741) we get for the correction, 



(vy = iv 5 



(3X 



2(2 



7T 




(3.75) 



where the function -F(£) is given by the integral, 



dX < - In 



A 2 -(l-0 2 



A 2 -( 



In 



A 2 -(l+e 



-1\2 



+ Jin 



1+0 2 

e + (a + 1) 



A 2 

2 



-e + (a - 1) 2 



In 



1-^-1)3 

-r 2 + (a + 1) 2 



<- 2 + (A - 1) 



(3.76) 
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which can be computed explicitly and is a combination (rather long) of log 
and poly-log functions. The integration over the main branch of the log is 
assumed. F(£) has the following properties: It is invariant with respect to 
inversion: £ — > Also F(0) = F(oo) = and it has a maximum at one. 
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4 Spin chain gas at low temperature 



Consider the Hamiltonian (13.41) in the regime when the inverse temperature 
(3 is large. In this regime the commutator term in the Hamiltonian becomes 
dominating over the quadratic part. If at the same time the Yang-Mills 
coupling is small enough, the one loop term becomes dominant also over the 
higher loop contributions. As discussed in [35] (see [36] for the background), 
at the large value of coupling in front of the commutator term causes configu- 
rations with non-vanishing commutator to be statistically suppressed. In the 
extremal case, when /3<? Y m ~~ * 00 > the ^ s f° rce d to remain in the valley of 
the potential i.e. the allowed space will be restricted to matrices satisfying, 

tr[f a ,f 6 ][* a ,* 6 ] = 0. (4.1) 

This condition is equivalent to restriction to commuting matrices, 

[*a,*d=0. (4.2) 

Since ty a are not Hermitian, the condition (14.21) is not enough for them 
to be diagonalizable. The Gauss law constraint, 

G=[tt a ,tf a ] = 0, a =1,2 (4.3) 

implies, however, that the commuting matrices are also normal. Even more 
generally, they should satisfy, 

[*a, = 0. (4.4) 

The condition ( 14.41) is already sufficient for simultaneous diagonalization of 
all \l/ a as well as \l/ a whose eigenvalues are complex conjugate to those of \I/ a . 

Passing to description in terms of eigenvalues of \fr a as well as ^ of 
ty a in quantum theory can be associated with change of variables in the path 
integral. The Hamiltonian takes the form, 

H = /itr^ a + \x ■ aa6tr* a * 6 = fi+ tr ^+^+ + (i_ tr (4.5) 

where in the last equality we used the eigenbasis of matrix (x ■ a). The 
change of variables to diagonal values and angular component gives rise to 
a Jacobian which is the product of the Vandermonde determinants for all \1/ 
and 

a 2 (^,vi/)= n n (4.6) 

a m>n 
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A common trick in the Schroedinger equation approach is to absorb this 
measure (see e.g. [37]) by a redefinition of the wave functional 



i> ^ yA 2 0,#)v>. (4.7) 

The continuity of the wave function requires that the sign of the root is 
chosen in such a way that after the redefinition the wave function becomes 
antisymmetric in diagonal modes. From the point of view of path integral ap- 
proach vanishing of the measure at such points means that the configurations 
with coinciding eigenvalues are prohibited since they have zero probability 
density. 

As large N or large L are concerned the number of excluded states is 
negligible with respect to the total number of states. Then, the partition 
function can be evaluated by the direct computation of the trace, 

Z(/i, x) = tr e ~ H ^ = e " M+ En N "~^~ E " N " = 

{7V±=0,1,2...} 
JV 

= exp (-JVpn(l - e -/1 +) - ln(l - e^-)]) . (4.8) 

As we discussed earlier, large L and S correspond to small values of the 
chemical potentials fi±. In this limit the partition function becomes, 

Z((j,, x) = exp (— iVln/i_|_//_) . (4.9) 

Solution to the Legendre transform equations (13.201) gives the following 
expression for the effective action, 

S eS (L, S) = N In L + L_ = iVln(L 2 - AS 2 ). (4.10) 



1 1 

1 - e"^+ 1 - e"^- 



5 String theory interpretation 

Let us give a qualitative string theory interpretation of above results. We con- 
sidered the effective action which is the thermodynamical potential counting 
the number of SYM operators of a given length L and given spin S. Ac- 
cording to the AdS/CFT dictionary this corresponds to the entropy of AdS 
(multi)string states having S unites of angular momentum in the (12) plane 
and L — S unites of the angular momentum in the (56) plane of 5* 5 , regarded 

14 Note, however, that in this case the redefinition is an operator action. 
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as embedded into M 6 . (This corresponds to the choice oc (0i + i0 2 ) and 
Z OC (05 + i0 6 )O 

Basically, we found three different regimes in the behavior of the model 
which are distinguished by different scaling of the thermodynamical quan- 
tities at large N: I. iV-independent, II. linear iV-scaling and, finally, III. 
quadratic iV-scaling. Recall that the meaning of N on the string side is 
the number of _D3-branes providing us with AdS 5 x S 5 background in their 
vicinity. Consider all three cases in detail. 

I. The effective number of degrees of freedom does not depend on N 
rather on L and S. This means that we are in a genuine stringy phase 
as strings with no internal structure present the effective fundamental 
excitations which we are counting. This phase occurs at moderate 
values of L: L< y/N. 

II. Linear iV-dependence says that we have a number of degrees of freedom 
proportional to the number of D3-branes. Therefore, it is natural to 
assume, that it is the D3-branes which are the fundamental excitations 
of the theory in this phase. This phase corresponds to large values of 
09ym an d moderate L. (To have the one-loop contribution leading one 
has to require also that (3g 4 is small.) 

III. The N 2 scaling, generally means that the Z)3-brane is not anymore 
a fundamental excitation but can be considered as a condensate or 
bound state of smaller objects. The situation is very similar to brane 
condensation in noncommutative gauge theories [38-41] (for a review 
see also [42]). This occurs at large L: L > y/N. Remarkably, the 
model in this regime is best described in terms of noncommutative 
gauge theory [19]. 

6 Discussion 

In this paper we considered a non-Hermitian matrix model describing the 
anomalous dimension spectrum of M = 4 SYM theory. In this study we did 
not recur to the integrability at infinite N neither to the description in terms 
of spin system. Moreover, N all the time was kept a finite statistically large 
number. We introduce a notion of temperature, which in the SYM theory 
plays the role of auxiliary parameter allowing restoration the density of the 
anomalous dimensions of composite SYM operators. At the same time due to 
the fact that the dilatation operator is the Hamiltonian corresponding to the 
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radial time in conformal theory this temperature corresponds to compactifi- 
cation of the SYM to the four- dimensional sphere 5* 4 . In the framework of 
string/ gauge duality this temperature is the genuine temperature of the dual 
string theory and can be attributed to the presence of a black hole [22]. 

For the two-matrix model we succeeded to compute the thermal partition 
function and the effective action depending on the total occupation number, 
which in the spin language is the aggregate spin chain length, and the occu- 
pation number difference, which corresponds to the total spin. We did this 
in the both approximation of high as well as of low temperature. The most 
attention and effort was given to the high temperature regime. In this regime 
we found a rather elegant approach based on the analytic expansion in terms 
of exponentials of chemical potentials. This approach allows an exact sad- 
dle point evaluation of the partition function. For sufficiently large chemical 
potentials fi± the contribution to thermodynamical potentials which scale 
like N 2 and N does not appear. This is compatible with the PET (Polya 
Enumeration Theorem) evaluations which holds for fi above the Hagedorn 
critical value. As in our case the spin is also included we may conjecture that 
the phase separation should occur along the critical line given b}0: 

1 - e-"+ - e - "- = 0, 

at which the first non-trivial mode to the saddle point equation appear. 
(In the case of zero x the above critical line reduces exactly to the Hagedorn 
criticality condition for /i, note that \i for the oscillator has the same meaning 
as the temperature.) In this regime the fundamental excitations are stringy- 
like. 

Using the incompressible liquid model approximation for the gauge field 
eigenvalue condensate we found a description of the model at small chemical 
potentials fj,±, which can be trusted at least qualitatively. For this regime we 
found a N 2 behavior, which corresponds to string bit phase. This phase is 
characterized by melting of strings into point objects: string bits. As a result 
the model looks as a system of N 2 interacting particles. 

At small temperatures the matrix fluctuations are bind to the diagonal. 
Therefore the effective number of degrees of freedom is proportional to N. As 
iV is the number of D3-branes in the string theory, the natural assumption 
is that they are the elementary excitations in this regime. 
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A Anti-holomorphic representation 



A very useful representation of wave functions of quantum oscillators is given 
by the anti-holomorphic representation. For the convenience of the Reader 
and to fix the notations we give here a brief account of this representation. 

Consider the oscillator algebra generated by the ladder operators a and 
a) which are subject to the commutation relations 

[0,0^ = 1. (A.l) 

We can represent the quantum states in the theory as anti-holomorphic 
function f(z) of a complex variable z. The action of the ladder operators on 
such states is given by 

a ] -»• z, d/dz, (A.2) 

i.e. the raising operator a) acts by multiplication by z, while the lowering a 
acts as a derivative. Obviously the commutation relation (lA.ip is satisfied 
by such definition. 

Now it is not very difficult to recoverer all formulas used to construct 
the representation. The oscillator vacuum is given by the zero mode of the 
lowering operator 

a |0) = = dipo/dz =>■ ipo — constant. (A. 3) 

We can choose (fo — 1. In this case the ortho normal set of the oscillator 
hamiltonian eigenstates is represented by 

\n) ~ 4=, (A.4) 



from which we can immediately extract the scalar product rule for arbitrary 
(anti-holomorphic) states described by f{z) and g{z): 

(f\9) = J ^e- Sz r(z)g(z). (A.5) 

A "well defined" operator F can be represented by an anti-holomorphic kernel 
F(z,z). The scalar product f]A.5|) implies the following rule for the kernel 
multiplication: 

F-G~ j e-™ w F(z, w)G(w, z) } (A.6) 
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In particular, the trace of operator F is given by 

tvF = J^e-^F(w,w)- (A.7) 

In many cases the normal symbol of an operator is used. Therefore, it is 
useful to know how to pass from one description to another. 
Consider an operator F represented in the normal form: 

F = F N (a\ a) = J2 fm,nW) m a n . (A.8) 

771,71 

Since the operators a) and a are normal ordered inside F^(a\a) one can 
treat it as usual function F^(z,z) of complex variable z. By analyzing the 
the matrix element of F between two oscillator states (m\ and \n) one finds, 
that the normal symbol Fn(ci,cl) is related to the anti-holomorphic symbol 
in the following simple way: 

F(z,z)=e 2z F N (z,z). (A.9) 

The last useful formula we want to give here is the anti-holomorphic 
symbol of the exponent of an operator with quadratic normal form: 

U = e~ pH , H = ha! a, (A.10) 

e.g. the anti-holomorphic symbol for the evolution operator of the harmonic 
oscillator. Again, the direct computation gives, 

U(z, z) = exp (ze- 0h z) . (A. 11) 

As a simple test let us compute the trace of ( 1A.11I) . Modulo vacuum 
contribution this gives the oscillator partition function, 

tre -PH = f ^ _ e _ ph , = 

J tt 1 — e~p 

The generalization to higher dimensions is straightforward. 

B Anti-Holomorphic representation with gauge 
symmetry 

Let us consider a situation when oscillator possesses a gauge symmetry, 

z i-f z 9 , geG. (B.l) 
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This pulls back as a G-transformation of the Hilbert space, 



^)H^) = r(9)^. (B.2) 
where T(g) is the generator of the induced action. Infinitesimally, 

T(g = e u ) « I + t(u), t(uM*) « ^(z e ") - ^), (B.3) 

where t(u) = t a (z)u a are the generators of infinitesimal gauge transforma- 
tions. Generally t a have the form, 

t a {z) = V a (z)8i, B i = -^i (B-4) 



where t l a (z) are just ordinary functions, 



du a dzj 



(B.5) 

9=1 



The gauging consists in projecting the Hilbert space to the subspace of 
invariant functions, which amounts to finding subspace 

t a -ip(z) = 0. (B.6) 

According to the canonical quantization procedure ( IB. 61) is equivalent to 
imposing the on-shell constraint condition since according to the second the- 
orem by E. Noether, operator t a correspond to first-class constraint in the 
classical theory. When G possesses an invariant Haar measure e.g. is a 
compact Lie group solution to flB.6j) is remarkably simplj^l 



^gauge invariant {z) = / dgif)(z 9 ), (B.7) 



G 

where we integrate over the gauge group using the invariant Haar measure dg: 
d{h~ l gh!) = dg. Hence the Haar integral ( IB. 71) can be used as the projector 
to gauge invariant subspace of the Hilbert space. It is not difficult to check 
that for an operator with gauge invariant kernel: K(z 9 , z 9 ) = K(x, z) action 
commutes with averaging over the gauge group, 

K - U-ip = I e-™ w K(z, w)^(w 9 ) = 



J dw^w_ e - WWR (- 9 ^ = Y[. K ■ if). (B.8) 
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I am grateful to Jeong-Hyuck Park for pointing my attention to this possibility. 
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In particular, the gauged analog of (IA. 12[) reads 



^ -rh f dwdw , r _ „ _ h -, 
trll • e p = / dg exp [-«j«j 9 + we w\ 

= Jdg [det(T^)-e- h )] _1 . (B.9) 

C The trigonometric sums 

We use the following trigonometric sums, 

mnk 

2^ cos (u9 m ) cos (u'8k) cos ((a; + a;') 9 n ) 

mnk 

— sin (u9 m ) sin (cu'Ok) cos ((a; + u/) n ) 

rrmfc 

+ sin (a>0 TO ) cos (uj'9k) sin ((a> + a;') 0„) 

mnk 

+ cos (c<j6> m ) sin (uj'Of.) sin ((a; + uo') 9 n ) 

mnk 

= Pu)Pu)'Pu)+u' — PujPuj'Plu+uj' + PujPuj'Puj+uj' + Pupu'Pu+u', (CI) 

where, 

= 2Jcosa;0 n , = ) y sin a;fl n . (C.2) 

n n 

As we are considering even distributions of 9 n , sums involving the sines 
vanish, i.e. p^ = for any u. Hence, the sum is reduce to, 

2j COS (uo9 mn + Uo'9 kn ) = PujPu'Puj+lu'- (C.3) 
mnfc 

D Incompressible liquid of eigenvalues 

Consider the static equilibrium problem for a system of particles with the 
pairwise interaction potential <p(X mn ) given by f)3.4ip where \ mn is the dis- 
tance between particles with numbers m and n. The structure of the potential 
is such that at small distance the interaction has repulsive character while 
at large separation it is attractive. Thus we may expect that N particles 
interacting with such a potential should have a static equilibrium position. 
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Let us further assume that this arrangement is compact for large N, i.e. the 
eigenvalue density vanishes exactly outside some finite region. Let us assume 
this region to be connected of size A. The condition of equilibrium requires 
that in the presence of non-zero density the potential $(A) defined as 

,A/2 

$(A)= / dvpirj^X-v), (D.l) 

J -A/2 

is constant $(A) = $o- Due to this condition the total energy is given by, 

E= fd\p(\)$(X) = iV$ , (D.2) 



where we used the normalisability condition for the density. 

Let us find the size of the distribution in the approximation of incompress- 
ible condensate. For this let us consider the Taylor expansion coefficients of 
the potential in the vicinity of origin: 

$(A) = $ = $(0) + $'(0)A + ±$"(0)A 2 + . . . (D.3) 

The zeroth term should give just $o ; 

/■A/2 

$o = p* / dri(p(ij), (D.4) 

J -A/2 

where is an average value of the density, in present approximation we can 
replace it by, 

P* = j- (D.5) 

The first term in the expansion vanishes automatically due to the sym- 
metry of the distribution. Consider the third term: 



,A/2 

d>£ = p, \ drjpirjMrj) = 2 P ^'(A/2). (D.6) 

J —A/2 



The vanishing of $q, which is zero pressure condition at the origin allows one 
to find the size of the condensate A, 

<pf(A/2) = 0. (D.7) 

Now, plugging the the potential (13.411) into into Eq. (ID. 71) we get: 



, N 

A = 2^fT, p» = (D.8) 

2^/jz^jn 
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i.e. the eigenvalue density is described by, 



N 



p(A)= J 2v ^T' V^- (DQ) 

[0, |A| > ^/Jl^JjZ. 

All this leads to the following result for the "energy" , 



M = -4^ (M^ (Jt) + Jf^ (yg) ^ (D,0) 



E Entropy of random spin states 

Let us compute the number of the states of spin s of a set of L spin 1/2 
states. It is instructive to compare the entropy of this set of random spin 
states with the one of the matrix model. 

The problem is related to one of computation of the multiplicity Vg of 
the irreducible representation of spin s in the expansion of the product of L 
spin- 1/2 representations: 



1/2 x 1/2 x ...1/2 = V^s, (E.l) 



L— times 



where we use the boldface letters s to denote the irreducible representations 
of spin s. 

To find the expansion (lE.lj) . let us observe that the product of an irre- 



ducible representation of s with 1/2 results in 

1/2 x s = (s+ 1/2) + (s- 1/2). (E.2) 

Consequent multiplication by 1/2 in the l.h.s. of spin 1/2 representation 
in (IE. II) can be represented by the diagram [El where each dot represent an 
irreducible representation of su(2). According to ( IE. 21) multiplication by an 
additional factor of 1/2 gives rise to a representation of representations of 
spin differing by ±1/2 with factor one each. On the picture this is represented 
by arrows. Thus the the total number of factors of spin s one gets after 
multiplication of L 1/2 factors is given by the number of distinct paths of 
length L one can reach the level corresponding to the spin s from the left-most 
position following the arrows. 

Thus, the problem is reduced to the computation of the number of dis- 
tinct paths to reach the level s. To do this let us note that the representations 
form families corresponding to straight lines starting from the bottom. Let 
us assign a number k — 0, . . . , L to each line. The multiplicity of irreducible 
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123456789 10 11 
L - 



Figure 1: Triangle of L products of 1/2 representations. The horizontal 
levels count (in the direction s) the spin of irreducible representations. The 
vertical dotted lined correspond to the length of the product. Thus the 
maximal length in this picture corresponds to the product of 11 1/2. The 
number to the left of each dot represents the number of distinct ways the dot 
can be reached from the "zero point" following the arrows. This number gives 
the degeneracy of the representation given by the level in the decomposition 
of the product of 1/2 in the number of factors is given by the horizontal 
coordinate L. 



representations found on a particular line satisfy simple properties. For ex- 
ample the multiplicities on the zeroes line are all equal to one. On the line 
number one the multiplicities are given by the sum of multiplicities from 
the zeroes line up to the next number of the sequence term. In general, the 
multiplicity on the fcth line is connected to ones on the (k — l)th one by the 
following recurrence relation 

n+l 

tf = iy*-i, "i = "2-i- ( E - 3 ) 

1=2 

The recurrence relation (IE.3I) can be solved which leads to the following 
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expression fo the multiplicities, 

(n + 2Jfe-l)! ,„ , 

k k\(n + k)\ K ' 

Now taking into account that k and n can be expressed in terms of the 
number of spins L and and the spin s as followJ^I 

L = n + 2k - 1, s = L/2 — k— (n— l)/2, (E.5) 

one can rewrite the multiplicities in the following form 

^' = ' 2 ' + 1 » (L/2-.)IW2 + . + l)! - (E ' 6) 

The check that the sum of the dimensions all irreducible representations 
taking into account the multiplicities is indeed 

]T(2 S + l)vV> = 2 L , (E.7) 

s 

is left to the reader. 

In the limit of large L one can use the Stirling's approximation to the 
factorials. This gives the log number of L random spin states of spin s to be, 

Random spin = L In L - \ {{L - 2s) hx{L - 2s) + (L + 2s) ln(L + 2s)}. (E.8) 

It is useful to compare this quantity to the spin entropy of the matrix oscil- 
lator in the main body of the paper. 
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